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We observe velocity-selective two-photon resonances in a cold atom cloud in the presence of 
a magnetic field. We use these resonances to demonstrate a simple magnetometer with sub-mG 
resolution. The technique is particularly useful for zeroing the magnetic field and does not require 
any additional laser frequencies than are already used for standard magneto-optical traps. We verify 
the efi^ects using Faraday rotation spectroscopy. 



PACS numbers: 42.50.Vk, 32.60.4-1 
I. INTRODUCTION 

Stimulated Raman transitions that couple atomic 
ground states with counterpropagating laser beams are 
resonant only within a narrow velocity band. This atomic 
velocity selection [T has proven to be a useful tool for a 
variety of experiments, including subrecoil Raman cool- 
ing [2], atom interferometry [5], and atom velocimetry [1]. 
Stray magnetic fields can adversely affect this process by 
shifting the magnetic sublevels, thereby perturbing the 
participating velocity bands [21 HI IS] . Conversely, when 
transitions occur between different magnetic sublevels of 
a single hyperfine level, velocity selectivity can provide 
an excellent measure of stray or applied magnetic fields. 

Elimination of stray fields to sub-milliGauss levels is 
particularly important for sub-recoil cooling processes [3J 
Typically, these fields are nulled by Helmholtz coils 
along each cartesian direction. Correct compensation 
currents can roughly be estimated by visual indicators 
such as atom expansion in an optical molasses, but 
these cues are strongly dependent on optical alignment. 
Stray fields can be directly measured using, for example, 
Faraday spectroscopy, which provides picoTesla sensitiv- 
ity [7j 8, 9j, but requires additional laser frequencies and 
time-resolved polarimetry. Measurement of vector mag- 
netic fields with magneto-resistive probes has been used 
for active compensation of both DC and AC fields, but 
needs several sensors placed externally to the vacuum 
chamber [5]. 

In this paper, we describe a simple imaging technique 
for measuring magnetic fields with sub-milliGauss reso- 
lution using a sample of cold atoms from a point trap. 
The technique relies on velocity-selective two-photon res- 
onances jl] (VSTPR) in a magnetic field, where the two- 
photon resonance occurs between different magnetic sub- 
levels within a single hyperfine level. When applied to 
atoms cooled in an alkali- vapor MOT, no additional laser 
frequencies are required, because the VSTPR pulse can 
be derived from the rcpumping laser beams. The two 
main requirements are 1) VSTPR beams along a hori- 
zontal axis and 2) a CCD camera whose optical axis is 
orthogonal to the propagation direction of the VSTPR 
beams. 




FIG. 1: Relevant energy levels for the Raman transitions. 
Magnetic sublevels are split by the Larmor precession fre- 
quency, lul- The polarization configuration is lin _L lin. The 
energy scale is exaggerated for clarity. 



II. BACKGROUND 

In this section, we briefly describe the process. Figure[l] 
shows the relevant energy levels. An atom moving with 
velocity v along the a;-axis is exposed to a light field 
composed of two counterpropagating laser beams with 
wave vectors ki and k2, where ki ~ —^2 — k along the 
X-axis. The polarizations of the two beams are lin _L lin. 
For an arbitrary B-field, this polarization configuration 
couples rrif — mi — Am = 0,±1,±2 magnetic sublevels of 
a single hyperfine level, where we choose our quantization 
axis along the magnetic field. Absorption of a photon 
from one beam and emission into the other results in 
a linear momentum change of ±/i(ki — k2) ~ ±2/ik = 
±2MwrX, where Vj. is the recoil velocity and M is the 
mass. The one-photon detuning, A is chosen to be much 
larger than the hyperfine splittings of the upper state. 

The two-photon detuning is defined here as (5 = cji — 
LJ2 — (Am)a'L, where lol is the Zeeman splitting. In a 
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FIG. 2: Experimental setup for observing VSTPR in a mag- 
netic field. 

small magnetic field, hiOL = gF^^BB, where g_F is the 
gyromagnetic ratio, and /is is the Bohr magneton. For 
s^Rb, gpf^B/fi = 466.74 kHz/Gauss [TU]. Two-photon 
resonance occurs for atoms whose velocity satisfies 

'5-^LS + ^D + 45,., (1) 

where the two-photon Doppler shift, Sd = 2k-v, and 
the recoil frequency Sr — hk'^/2M. The relative light 
shift, (5is is a weak function of the participating magnetic 
sublevels, which are both in the same hyperfine ground 
state. Apart from Sls, and with lui — uj2, the resonance 
condition is satisfied for atoms with v = vq ± such 
that 2k- Vo — {Am)ujL. 

Throughout the duration of the VSTPR pulse, reso- 
nant atoms oscillate between the two momentum states 
separated by 2/ik. Because the atoms are initially con- 
fined in a point trap, an image of the atom cloud after 
expansion is a spatial map of the average velocity distri- 
bution, which has been perturbed by the VSTPR pulse. 
A freely expanding cloud has approximately a smooth 
Gaussian velocity spectrum along x (the VSTPR beam 
direction), but the momentum oscillations that occur for 
the resonant atoms alter this average velocity distribu- 
tion. Images taken along a camera direction orthogonal 
to X record these narrow perturbations. 

III. EXPERIMENT SETUP 

The layout of our apparatus is shown in Fig. [2j The 
experiment begins with a vapor cell MOT containing 10^ 
*^Rb atoms. The MOT diameter is w500 /im and the 
temperature is «200 /iK. The cooling beams (detuned 
-18 MHz from the F—3 - F'—i transition) are derived 
from an extended cavity diode laser (New Focus Vortex, 
model 6013) that seeds a 120 mW laser diode (Sharp 
GH0781JA2C). 60 mW is directed to one port of a 2 x 
3 polarization maintaining (PM) fiber splitter. Each of 



the three output fibers carries 12 mW. The outputs are 
coUimated using 100-mm focal length, 50mm-diameter 
achromats, giving 1/e^ beam diameters of 24 mm. The 
three beams propagate along orthogonal directions and 
are retrorefiected; one pair is vertical and the other two 
are in a horizontal (x-y) plane. The repump light, con- 
necting = 2 ^ i^' = 3, is derived from an independent 
Vortex laser that also seeds a diode. In normal MOT op- 
eration, 15 mW of repump light is coupled into the other 
port of the 2x3 coupler. 

Our VSTPR beam is spatially filtered by PM fiber, 
and is coUimated by a 60mm gradient-index singlet lens 
(1/e^ beam waist ujq — 7.5 mm). We use up to 20 mW 
laser power. It is retrorefiected in a lin ± lin configura- 
tion. For B-field control, we use three orthogonal pairs 
of Helmholtz coils. The magnetic field at the atom cloud 
has components Bi = OLi{Ii — iQi) where a.^ are the slopes 
dBi/dIi, li are the applied currents, and I^i are the cur- 
rents required for compensation along each Cartesian di- 
rection. The VSTPR beam travels horizontally along 
the axis of the x-directed coil pair. For all results in 
this paper, we have used a second beam derived from the 
repump laser as our VSTPR beam (A = 3GIIz). 

At time T = 0, the atoms are released from the MOT 
by extinguishing all laser beams and the MOT coils. The 
bias magnetic coils remain on. We do not perform any 
molasses cooling, because the large velocity spread of the 
hotter sample of atoms provides greater range over which 
velocity selection can occur. At time Tr ~ 15 ms, the 
VSTPR pulse is switched on for 5 ms and at Ti = 40 ms, 
the MOT cooling and repump beams are switched on to 
image the expanded cloud onto the CCD camera. 



IV. RESULTS 

We control the magnetic field by changing the current 
in the Helmholtz bias coils. We first demonstrate the ef- 
fect by changing the current in the ^-directed Helmholtz 
coils. Fig. [3|i shows typical images recorded by the 
CCD camera for four different current settings. The 
VSTPR pulse creates perturbations that appear as verti- 
cal stripes in the expanded cloud. Also shown in Fig. [3|l 
are the cross sections after subtracting the images ob- 
tained with no VSTPR pulse. The spatial location of 
each peak corresponds to the average velocity class satis- 
fying 2k-vo — ±(ATO)a;i that participated in the VSTPR 
process. In general, lul is a function of position due to 
spatially varying magnetic fields. These gradients could 
tip, bend, or blur the resonant stripe features. Under 
normal operating conditions, however, we do not observe 
these effects. With the B-field primarily along the z- 
axis and using lin _L lin polarizations, the VSTPR pulse 
mainly connects Am — 1 transitions. The cross-sectional 
profile of each vertical stripe can be estimated by several 
functional forms. For simplicity, we fit each peak to two 
Gaussian profiles such that the total area is 0, but we 
describe a more exact fitting function below. 
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FIG. 3: a) Left: Raw images of the expanded atom cloud 
at different Helmholtz currents after exposure to the VSTPR 
beams. Right: Corresponding cross-sectional profiles after 
background subtraction. Fits are shown as solid lines, b) 
Larmor precession frequency extracted from the stripe sep- 
aration as a function of current in the z-directed coils for a 
nonzero By field (circles) and for Bx = By = (diamonds). 
Fits shown as solid lines. 



In addition to the peaks corresponding to velocity 
classes at ±vo, there is a peak at v = due to Am = 
transitions that arises from the longitudinal magnetic 
field component. This is a useful marker for balancing 
the overall velocity distribution. We rarely observe fea- 
tures corresponding to Am = ±2 transitions. Detailed 
calculations of the relative transition strengths will be 
published elsewhere, but in general, the strengths of the 
Am = ±2 transitions are roughly two orders of magni- 
tude smaller than those for Am = 0, ±1 for A = 3 GHz. 
It is important to note that the appearance of narrow fea- 
tures in the expanded cloud is not the result of cooling, 
because there is no dissipative force. For smaller A, on 
the order of the hyperfine splittings of the D2 manifold, 
spontaneous scattering events can lead to magnetically- 
induced laser cooling [11]. 



The z-directed Helmholtz bias coils are each made of 
24 gauge wire wound on an 8" vacuum flange, 1" long, 
and centered 5.8 cm from the MOT, producing a field 
of approximately 1.5 G/A. Our fits to Fig. [3] show that 
after 35 ms of falling time, the separations of the stripes 
for the 4 current settings 260 mA, 330 mA, 400 mA, 
and 450 niA are 3182(4) /im, 4538(4), /im, 5902(4) ^m, 
and 6873(4) /im. The listed errors are statistical. Our 
leading systematic source of multiplicative error is the 
pixel calibration on the camera of «0.2%. In section VI 



we briefly describe a second calibration procedure that 
removes pixel calibration errors. The 4 /im error in stripe 
separation corresponds to an error of 300 /iG. Note that 
the error in magnetic field will be explicitly dependent on 
the atom species through gp and the VSTPR wavelength. 

The scalar magnetic field measurements for several cur- 
rent settings are shown in Fig. |3]d. Because the stripe 
separation is proportional to {B^+ By + B^)^^^, a plot 
of ujl versus the current in the z-directed bias coil traces 
a hyperbola, the minimum of which determines the field 
component perpendicular to z and the compensation cur- 
rent Iqz- We show two cases, one with nonzero By, and 
one for the case in which B^, By have been zeroed. From 
a fit to these plots, we extract = 1.524(2) G/A and 
Iqz — 243.1(2) mA, which corresponds to a compensation 
level of 300 fxG. In practice, it is simple to zero the mag- 
netic field by viewing real-time images of the expanded 
cloud and adjusting the currents along each axis for min- 
imum stripe separation. We note that when the total 
magnetic field is close to zero, the stripes begin to over- 
lap and are no longer resolved. Compensation is achieved 
when the overlap is maximized, resulting in a single nar- 
row feature. In our experience, this real-time adjustment 
of the stripe separation results in compensation to milli- 
Gauss levels without any data analysis. 

The visibility of the stripe features is dependent on a 
few factors. First, since the image on the CCD camera 
is a convolution of the initial MOT size with the velocity 
distribution, the contrast increases for trapped samples 
with smaller physical dimensions. Optimally, the imag- 
ing should be performed after the cloud has expanded 
enough that two velocity classes separated by 2/ik/M can 
be resolved. If the initial MOT has a radius R, this means 
that the imaging should be performed after a time R/vrec 
from the release of the atoms from the trap. In practice, 
the features are easily observed with imaging times sig- 
nificantly less because the effect does not require that 
the recoil velocities be resolved, only that perturbations 
to the average velocity distribution can be observed. 

Other parameters that control the visibility of the 
stripes are the duration and timing of the VSTPR pulse. 
Because this measurement is time-averaged over the du- 
ration of the pulse, shorter pulses reduce blurring effects 
due to time- varying fields. Furthermore, if applied at 
— Til 2 they can maximize stripe contrast. Figure |4|l 
shows the effect on stripe contrast as a function of 
for atoms in zero field. An atom initially at the origin 
will return to the origin if its momentum is reversed by 
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FIG. 4: a) Importance of timing of the VSTPR pulse. An 
atom moving with hk has its momentum reversed at Tr- If 
this pulse occurs at Ti/2, atoms initially having iLhk are over- 
lapped at Ti. b) If the pulse arrives later or earlier, contrast 
is reduced, c) Pictures of a single stripe using different Tr, 
showing optimum contrast at Ti/2. d) Stripe cross-sections 
for Tr — 0.5Ti and T,. — 0.25Ti, along with fits using Equa- 
tion [2] (solid line) . 



a TT-pulse at T,/2. For AT = - TJ2 ^ 0, the atoms 
are still deflected, but the atoms with initial momenta 
of ±hk no longer spatially overlap at Ti (Fig. |4|d). In 
Fig. |4j:, we show this effect experimentally using VSTPR 
pulse durations of 200 /zsec at different T^. This simple 
geometrical picture suggests an approximate functional 
form for the background-subtracted stripe cross section 
under these conditions: 



Y{x) = G{xQ-2vrAT) + G{xo + 2vrAT) 

- G{xa - VrTi) - G{xo + VrTi) (2) 

where G{xo) is a Gaussian centered at xq- In Fig.jiji, we 
show fits when Tr = O.bTi and Tr = 0.25Ti. For all other 
data presented in this manuscript, we used pulse dura- 
tions of 5 msec. Although AC magnetic fields were not 
compensated, these longer pulses showed no measurable 
broadening in our experiment. The 7r-pulse duration de- 
pends on the particular magnetic sublevels involved, so 
for our spin-unpolarized sample we generally choose a 
pulse duration that provides consistently strong signals 
over a broad range of Raman pulse intensities. For a typ- 
ical VSTPR beam intensity of 10 mW/cm^, the resonant 
2-photon Rabi frequency is «27rxl0 kHz. 
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FIG. 5: a) Schematic of Faraday spectroscopy setup. WP: 
WoUaston prism, b) Typical Faraday signal, c) Compari- 
son of magnetic field measurement with VSTPR (circles) and 
with Faraday spectroscopy (squares) for similar conditions as 
Fig. |3]d. Solid line is a fit to the Faraday data. Dotted line is 
a fit to the VSTPR data. Dashed line indicates the minimum 
of the hyperbola as measured by Faraday spectroscopy. 



V. COMPARISON WITH FARADAY 
SPECTROSCOPY 

We have verified this VSTPR technique by using 
Faraday rotation spectroscopy to measure the magnetic 
field f7J [SI H]. To perform these measurements, an ad- 
ditional pair of laser beams is used along the cc-axis 
(Fig. [5^). The atoms are optically pumped into the 
F — 3, nip = 3 stretched state (in the x-basis) by a 
100 /iS (T+ pulse connecting F — 3 ^ F' = 3. This beam 
contains a small amount of repumper light to keep the 
atoms in F=3. When this light is extinguished, the atoms 
begin precessing freely. A linearly polarized probe beam 
with ~100 fiW and waist ujq = 500/im passes through the 
atom cloud to a simple polarimeter consisting of a Wol- 
laston prism that splits the probe beam into two orthog- 
onal polarization states that are detected by a balanced 
photodetector [12 . We make these measurements at the 
same time delay as the VSTPR measurements. A typical 
Faraday signal is shown in Fig. [5|d. The comparison of 
the VSTPR and Faraday techniques is shown in Fig. [sj;. 
Because we cannot compare Faraday and stripe measure- 
ments near zero field where the stripes are unresolved, we 
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FIG. 6: a) Image of atom cloud (after background subtrac- 
tion) taken with several sidebands separated by 100 kHz im- 
posed on one VSTPR beam, b) Cross sectional profile. 

show in Fig. [5|: the stripe data with and without a trans- 
verse field along y. This transverse field allows a stripe 
measurement at Iqz, which is in good agreement with 
the Faraday measurement (dashed line in Fig.jsj:). From 
the Faraday measurements with no transverse field, also 
shown in this figure, we derive = 1.547(3) G/A and 
Iqz = 242.5(1) mA, both of which agree well with the 
VSTPR technique. 

VI. CALIBRATION 

To obtain correct values of the magnetic field, the 
stripe separation must be carefully measured. Without 
calibration, the technique is still useful for determining 
the loi of each coil by simply finding the minimum stripe 
separation, which is independent of this source of sys- 
tematic error. Because this is an imaging technique, 
good estimates of spatial calibration can be made sim- 
ply by measuring the magnification on the CCD camera 
if high-quality imaging lenses are used. Most lenses ex- 
hibit some degree of aberrations that make accurate pixel 
calibration difficult below the 0.5% level. Even without 
this error, other slight systematic errors, such as the ex- 
act functional form used to fit the stripe cross section 
may remain. In this section, we describe a technique for 
calibrating the Zecman shifts indicated by the stripes in 



a more direct manner that eliminates these systematic 
errors. 

Instead of relying on an accurate pixel calibration, the 
splittings of the stripes can be determined by making 
uji = u)2 + Si2, where S12 is a frequency shift imposed 
by an RF source. The resonant velocity classes are now 
determined by 2k- v — (Am)^^ + ^12- To demonstrate 
this idea, the retroreflecting mirror in Fig. |2]is replaced 
by a counterpropagating beam of the same diameter and 
power. This counterpropagating beam is derived from 
the original, so it is phase locked to loi , but its frequency 
is shifted by two acousto-optic modulators (AOM) to 
achieve small 5i2 (<1 MHz). Additionally, we frequency 
modulate one AOM so that its instantaneous frequency 
is ionp + Asin^LUmt), where uhf is the drive frequency, 
ujm is the modulation frequency, and A is the maximum 
frequency deviation. This imparts frequency sidebands 
at ncurni where n is an integer, so that multiple S12 are 
produced simultaneously. In Fig. |6] we show an image 
and cross section taken with uim = 27r x 100 kHz. For 
612 = 0, the range of measurable Zeeman shifts is limited 
to the Doppler width of the atom cloud (~1G). A nonzero 
S12 overcomes this limitation by shifting the stripe to an 
accessible velocity class. 



VII. CONCLUSION 

We have used velocity-selective resonances between 
magnetic sublevels of a single hyperfine level in ^^Rb 
to measure magnetic fields in a cold atom cloud. The 
resonances are easily observed with no additional laser 
frequencies than are required for MOTs, and can be used 
to measure magnetic fields with sub-mG resolution. Be- 
cause of its simplicity, this technique should prove espe- 
cially useful for aiding magnetic field compensation, for 
which purpose no calibration is required. 

This work was funded by the Defense Advanced Re- 
search Projects Agency and the Office of Naval Research. 
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